(N 



New Applications of Quantum Algebraically Integrable Systems in 

Fluid Dynamics 

Anne Boutet de Monvel, Igor Loutsenko, Oksana Yermolayeva 

Institute de Mathematiques de Jussieu, UPD Paris VII 
Institute Henry Poincare, UPMC Paris VI 

e-mail: aboutet@math.jussieu.fr, loutsenko@math.jussieu.fr, yermolay@ihp.jussieu.fr 



O '■ Abstract 

Quantum algebraically integrable systems are non-trivial generalizations of Laplacian 
operators to the case of elliptic operators with variable coefficients. We study corre- 



sponding extensions of Laplacian growth connected with algebraically integrable sys- 
tems, describing viscous free-boundary flows in non-homogenous media. We introduce 
a class of planar flows related with application of Adler-Moser polynomials and also 
construct solutions for higher- dimensional cases, where the conformal mapping tech- 
nique is unavailable. 

1 Introduction 

We start with an inverse problem in fluid dynamics. Before posing the problem, we explain general 
settings and notations. 

Consider a viscous flow or a flow of liquid in a thin (not necessarily planar) layer of non- 
homogeneous porous medium. The layer can be viewed as a 2-dimensional surface embedded 
in the 3-dimensional Euclidean space. We let the layer curvature, permeability, porosity and 
thickness depend on the surface spatial coordinates x, y. We can choose x, y such that locally 



dl 2 = G(x,y)(dx 2 + dy 2 ) 

where dl is the surface length element. The surface area element is <i£ = Gdxdy, and the volume 



of the liquid that can be absorbed in the range x + dx,y + dy equals 

dV = rjhdT, = r/hGdxdy, 



where rj = t](x, y),h = h(x, y) are the medium porosity and the layer thickness, respectively. 

In the porous medium, the flow velocity v = (dx/dt,dy/dt) is proportional to the gradient 
V = (d/dx, d/dy) of the pressure P 



v = ~VP, 

Vg 



where y) is the medium permeability. 
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It is seen from the above that only the two combinations of variable coefficients, namely 



r]hG, —= 

Vg 



enter the flow equation of motion, and it is convenient to absorb h and G into definitions of the 
other coefficients. Therefore, without loss of generality, we can consider the flow in the plane 
parametrized by the complex coordinates z = x + iy,z = x — iy, choosing rj and k to depend 
on z, z, while setting remaining coefficients to unity It is also convenient to redefine pressure 
P-t-P. 

Then the liquid volume conservation leads to the continuity equation 

(V -r)v) = 0, (1) 

while the dynamical law of motion rewrites as 

v = kS/P. (2) 

We consider a situation where the liquid occupies a bounded, simply-connected open region Q of 
the plane, whose time evolution fl = fl(t) is induced by the flow. 
At fixed time t, the pressure is constant along the boundary 

p(dsi(t)) = p Q (t). (3) 

Note, that the in the case of the simply-connected domains considered here, the dynamics is 
independent of Po(t), and for convenience the latter can be set to zero. 

The normal velocity of the boundary v n and that of the flow coincide at dfl 

v n = n ■ v if z G cftl (4) 

The flow is singularity driven. For instance 

Zl)Vv z li z l) ~[ \\ z ~ z l) \ z ~ z l) J 

when k + 1 multipole sources are located in vicinity of z = z\ in interior of the fluid domain Q. 

Equations flTJ - (J5]) set a free boundary problem where evolution of the boundary dQ(t) is com- 
pletely determined by the initial condition dQ(0) and dynamics (strengths) q(t) = q(t), Hj{t),j = 
1, .., k of the sources. 

It follows from equations ([T]), (j5J) and (JSJ) that the pressure satisfies the following linear PDE 

(VktiV)P = 2irq[8(x - x x )8{y - y x )) (6) 
where q is a differential operator of the order k 



j=i \ ' 
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and 5 denotes Dirac delta-function. 

Since permeability and porosity are always positive, the operator in the lhs of the last equation 
is elliptic and this class of free boundary problems is called elliptic growth. For the case of 
homogenous media 

K = const, i] = const, (8) 

the above operator is Laplacian and a subclass of related growth processes is called Laplacian 
growth or Hele-Shaw problem (see e.g. [11]) . 

Inverse problem: Consider an inhomogeneous medium with given permeability k(x) and 
porosity rj(x). In this problem one wishes to control the boundary dynamics using finite number 
of closely placed fixed sources/sinks by varying their strengths in time. 

The solution to this problem may be provided by finding a number of sources (or sinks) 
that better approximate the dynamics in quest. The shape of domain is approximated by time 
dependent polynomial conformal mapping 

i 

z = F(w,t) = z 1 + r(t)w + s ^u i (t)w\ f = r, (9) 

i=2 

from the unit circle \w\ = 1 in the "mathematical" plane w. For the case both permeability 
and porosity being functions of one space variable the problem can be solved using the following 
number of sources 

k + 1 = (deg(3(x)) + n(n 2 + ^ + deg(FH). (10) 

Here n is a number of Adler-Moser polynomial p n (x) and S(x) is a polynomials with free coefficients 
that approximate characteristics of the medium inside fl: 

KV = £(x)- 2 , £(x) = Pn X , v 7]{x) =p n - 1 (x)S(x) (11) 

Pn-l{X) 

The intrinsic property of Adler-Moser polynomials p n (x) consist in possessing n + 1 real free pa- 
rameters, two of them coming from shift and scaling symmetry and other n—1 defining polynomial 
coefficients. This, together with freedom in definition of S(x), allows approximation the medium 
characteristics. Increase in the quality of approximation will, in general, result in increase of the 
number of sources as seen from ( ITUI) . 

Example: Let us consider simplest example with linear porosity rj(x) = x. In the case the 
first two Adler-Moser polynomials 

Pi = x, p 2 = x 3 + ti 
can approximate with k(x) close to the following, 

x 

K= (x 3 + tx) 2 
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Then, according to (|T0l) . 

k + 1 = 4deg(z(w)). 

For instance, 4 sources are needed to sustain a circular form of domain, i.e. z = Z\ + r(t)w, for 
k(x) of the present example. 

The above problem and examples served as motivation for the present paper that have emerged 
in relation with our recent studies in quantum algebraically integrable systems (QAIS) and in 
elliptic growth. 

In the sequel we will establish important connections between the theory of QAIS, integrable 
elliptic growth and generalized quadrature identities. Later, we will consider higher-dimension 
cases, when the conformal mappings are inapplicable and the result is uniquely due to advances 
in the theory of QAIS. 

2 Quantum algebraically integrable systems as generalisa- 
tions of Laplace operators 

Let us review some common facts on QAIS. The classical Liouville integrability is equivalent to 
existence in a Hamiltonian system a maximal set of Poisson commuting invariants: For a given 
classical integrable system with n degrees of freedom there cannot exist more than n functionally 
independent integrals in involution. 

In the case of quantum integrable system in n dimensions, integrals form a commutative rings 
of differential operators generated by n independent elements and such a ring is called complete. 
Algebraic quantum integrability exhibits when there appears an additional, n + 1th independent 
generator of the ring of commuting differential operators. Such rings of commuting differential 
operators are called overcomplete. 

The simplest example of nontrivial QAIS comes from the theory of Calogero-Moser operators 
in one dimension. This example is a quantum system with the Hamiltonian 

h = dl - - 2 . 

x z 

For our purposes, it is convenient consider a gauge-equivalent differential operator 

1 2 

L\ = xHx~ l = x 2 d x —d x = d 2 x d x 

x A x 

For this system, there exists another differential operator 

M x = d x - - 2 d x + \d x 

that commutes with Li and does not belong to the ring generated by L\. Thus, in one dimension, 
we have the overcomplete commutative ring generated by two independent elements L\ and M\. 
It is important to note that this ring is not trivial in a sense that it is not a subring of any ring 



4 



of commuting differential operators generated by a single element, since there is no a first-order 
integral for the above system. 

The following example concerns the 2-dimensional case. In d = 2 the ring may be formed by 
the following set of commuting differential operators: 

Li = x 2 V^V = x 2 d x ^d x + d 2 , M 1 = dl i -^d x + ^d x , K 1= d 2 . (12) 

Let us note, that operator L\ is not trivial in the sense that it is not equivalent, up to change 
of dependent and independent variables and multiplication, to any Beltrami-Laplace operator on 
some 2-dimensional manifold. The operator Li describes elliptic growth with ktj = 1/x 2 in ([6]). 

Another important concept of algebraic integrability is the intertwining operator. Consider 
the differential identity 

AT = TL (13) 

that relates the Laplace operator with a differential operator of a second order L. Operator T is 
called the intertwining operator and (JTBl is an intertwining identity. 

The simplest example of intertwining identity refers to factorization of the one-dimensional 
Laplace operator d 2 : 

d 2 ( 1 d \ ( d 

x— 1 



dx 2 \ x dx J \ dx 
and 

d „ \ „, / d \ ( 1 d \ ( d 



X dx y dx y \xdx) \ 'dx 

from which we obtain the simplest intertwining identity 



7\A = L X T X T 1 = x^-1 (14) 



that relates the Laplacian with the Calogero-Moser operator L\ of ( fl2j) . 
Higher intertwining identities 



T A — IT 

tl — ±J n ± n 



involve operators of higher order, like the following examples: 

, , d n\ ( d n-l\ ( d 1\ T 2 „_ 1 „ 

T n = x n ( ] ( )...[- ) , L n = x 2n V^V 



OX X J \ OX X I \OX X 1 X 



An 



The main purpose of the intertwiner is to produce eigenfunctions of the interlaced operators from 
eigenfunctions of Laplacian: 

A0 = A0 LT[<j>\ = \T[<j)} 

Classification of QAIS is an open problem. Presently, algebraic integrability is found in systems 
related with finite- reflection Coxeter groups and their special deformations [TOj [7J [6]. In two 
dimensions, most general class of QIAS is related to soliton solutions of the KDV hierarchy [THE], 
while the algebraically integrable systems with coefficients depending on one variable only are 
related to rational solutions of the KDV hierarchy (Adler-Moser polynomials) [T3~l IT2] . 

In the following section we review applications of QAIS to generalization of harmonic analysis 
emerging in free-boundary problems of fluid dynamics. 



5 



3 Mean value integral identities, conservation laws and 
algebraic domains 

Let <j)(z, z) be a time-independent function satisfying 

Vnr]V<p = 0, zett (15) 

in whole Q, including z — Z\. 

Let us now estimate the time derivatives of the following quantities 

M[(f>] — / r](pdxdy. 
Jn(t) 

Considering an infinitesimal variation of the fluid domain Q(t) — > Q(t + dt), we get 

dM[cj>} r 

— — = <p v n r](pdl, 
dt Jdn(t) 

where dl is the boundary arc length. From ([2]), (E]), fll]) it follows 

^MM. = I (P Kr] V<f) - <j)KT]VP) ■ ndl. 
dt Jdn(t) 

Applying the Stokes theorem and remembering that P and <j) satisfy (??), ( |T5l) . we get 

dM[(j)} 



dt 



27rf MCsx.Zi). (16) 



Note that mixed derivatives are absent in q (c.f. ((?])), for by f|T5|) . is expressed through first 
derivatives of 4>. cjj is the complex conjugate of qj, since both <p(z, z) and (fi(z, z) satisfy f fT5|) . 

It follows that M[(j>] is conserved for any solution of f[T5"|) . a such that q*[(j)}(zi, Z\) = 0. 

The conservation laws have been first obtained for the homogeneous medium flows in the 
variable-coefficient generalization seems to be first presented in [8]. 

The Laplacian growth ([8]) is the simplest example, where the conservation laws can be written 
down explicitly [9]. In this example, any (anti) analytic in Q function satisfies (|T5|) 



<t>{z,z) = f(z) +g(z), for k = 1, rj = 1, (17) 
where f,g are (anti) analytic in Q and the quantities 

(f(z)(z - Zl ) k+1 + g(z)(z - Zl ) k+1 ) dxdy 



n(t) 

are integrals of motion for the free-boundary flows driven by a multipole source of order k located 
at z — Z\ in homogeneous medium. 



6 



Returning to the general case, we integrate ( IT6l) getting 

M[ ( f>)(t) = M[<f>)(0) + 27rQ[(f>}(z 1 ,z 1 ), 

where 

r* f d j d j 

Q = j a ?(W ' = « + E(<3iw + *aF 

Therefore M[0](t), and consequently a form of the domain, does not depend on the history of the 
sources and is a function of "multipole fluxes" 

Q = [ q(t')dt', Qj = [ qj (t')dt', Qj = [ qj(t')dt', j = l..k 
Jo Jo Jo 

injected by time t. 

Now consider the special case when M[0](O) = that describes the injection of the fluid to an 
initially empty medium. In such 

r](z, z)(j)(z, z)dxdy = 2rcQ[(f)](z 1 , zi), (19) 

Equation (fl9|) is a generalization of mean value or quadrature identities appearing in the theory 
of harmonic functions [TTj to the case of elliptic equations with variable coefficients. The simplest 
example of a quadrature identity is a mean value theorem for harmonic functions. Special do- 
mains for which the quadrature identities hold are called quadrature domains in the theory of the 
harmonic functions. 

Note, that above derivation can be also adapted to higher dimensional case where similar 
identities hold. 

Now, we return to our inverse free-boundary problem stated in the introduction. It translates 
to the language of quadrature domains as follows: Given a domain approximated by polynomial 
conformal mapping (Q, find k and r\ that best approximate given permeability and porosity and 
for which a quadrature identity holds in the above domain. 

Below, we are going to consider the case of stratified media, when permeability and porosity 
depend on one coordinate. 




4 Planar free-boundary flows in stratified media and ra- 
tional solutions of the KDV hierarchy 

Below we consider solution to the inverse problem for stratified medium presented in the intro- 
duction. 

Theorem 1 Let Q(t) be a domain in the plane given by a polynomial conformal mapping of the 
unit disc (TJ|). Then this domain can be formed by injection of fluxes Q and Qj,j = l..k through a 
multipole source of order k at z = z\ in a medium with permeability and porosity given by 07]) . 
The order of the source is given by l[TU\) . 
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The proof is based on a constructive procedure of derivation of a linear system of equations for 
values of fluxes Q, Qj that is a verification of quadrature identity ( 1191) for all regular in Q solutions 
of the elliptic equation 

e(x)2V ^ V0 = ° (20) 

with £ given by ( ITTj) in an algebraic domain defined by mapping (j2J). The "x- dependent" part of 
the above elliptic operator £(x) 2 d x £(x)~ 2 d x is gauge equivalent to a Schrodinger (Lax) operator 
connected with rational solutions of the KDV hierarchy with potential expressed through Adler- 
Moser polynomials [12]. 

The proof follows general argumentation that was used for system with finite reflection invari- 
ance in [I]: Here we just outline a procedure needed for finding values of fluxes. 

In the case of equation fl20l) . the intertwiner between the elliptic operator in (120]) and Laplacian 
(cf (|13p ) can be expressed in a form of a Wronskian 



T[f]{Xiy) =mMx),---,Mx)J(x,y)] 

where the Adler-Moser polynomials p n are also expressed as Wronskians 

p n (x) = Wliptix), . . . , ip n {x)) 

of functions 

ip i+ i(x)" = i/ii(x), ^i(x) = X 

The nth Adler-Moser polynomial depends non-trivially on n — 1 free parameters ("KDV times") 
that emerge as integration constants for ipith in the above equation. 
It then follows that general solution of ( 12"U|) can be written down as 

cf>{x,y)=T[f(z)+g(z)} 

with T given by (12T1) and f(z),g(z) being functions (anti) analytic in Q. Taking ffTTl) into account 
we get 

# = S(x)W[Mx), . . . , $ n {x), f(z) + g(z)\ 

which is linear differential operator with polynomial in z,z coefficients acting on f(z) or g(z). It 
is sufficient to consider its action on f(z): 



v4> = ^2Pj(z,z) 



dzi 



i=0 

We can now evaluate lhs of the quadrature identity ffl9l) . Using Green theorem, we get 



/ 



ncpdxdy = I Rj(z, S)^f-dz, ^p^- = ^Pj(z, z) (22) 

j—Q J dQ> 
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where 

R j = c ^ 1 + de s(5)+i0+ 1 )/2 + l OW e r order terms in z, z (23) 

is a polynomial of deg(S') + 1 + j(j + l)/2 th power in z, z. 

Since dQ is an image of unit circle under the mapping (jHJ), and since on the unit circle 

w = l/w, F(w) = F(l/w) 

we can rewrite (1221) as integrals taken around the unit circle 

f V cf>dxdy = V / F'(w)R j (F(w),F(l/w)) (-^d w Y [f{F{w))]dw. 
n ,_ „ J\w\=i \f \ w ) J 



j=0 



The map (Q is analytic and polynomial with F'{w) / on the unit disc, and since Rj(z,z) are 
polynomials in z, z, the above integrals can be evaluated as a sum of a finite number of residues 
at w — 0. Taking fl23|) and ([9]) into account, we get 

(deg(S)+l+n(ra+l)/2) dog F+n-1 



w=0 



Now we consider rhs of quadrature identity fll9p . According to (ITS]) we have 

Q[<P](z 1 ,z 1 ) = Q[T[f(z)]] z=Zl = 

j=l v 7 w =0 j=l m=0 v 7 

Equating the above expressions for lhs and rhs of the quadrature identity, we get k+n + 1 equations 
for 2k + 1 unknowns Q, Qj, Qj,j = l..k with k given by ( FlOl) . The condition that Qj is a complex 
conjugate of Qj and Q is real uniquely defines all fluxes. 

Concluding consideration of highly viscous planar flows, it is interesting to note that Adler- 
Polynomials, and rational Baker- Akhieser function of the KDV hierarchy in general, seem to be 
ubiquitous in integrable problems of two dimensional fluid mechanics: They also emerge on the 
other extreme, in the theory of inviscid flows that are finite-dimensional reductions of the Euler 
equations, describing equilibria and uniform motion of vortices in two dimensions [T^fl [T5| IT6] . 



w=0 



5 Higher-Dimensional Flows 

Three-dimensional viscous free-boundary flows are also of practical interest, while in dimensions 
four and higher the interest is mainly due to generalization of theory of harmonic functions and 
mean-value quadrature identities. 

Since the conformal mapping technique is unavailable here, we propose to write-down multidi- 
mensional quadrature identities using methods of construction of fundamental solutions for QAIS 
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Let us demonstrate application of these methods by a direct approach using example of an 
expanding sphere (Ci - C() 2 + (C2 - Q 2 + ■ ■ ■ + (Cd - Q 2 = r{tf centered at C' = (& • • • , Q in 
d- dimensional Euclidian space. The flow is governed by the following elliptic equation 

C[P](C) = tr^mCi - CO • ..S(Cd - C)L (24) 
where <Jd-\ is area of d — 1 dimensional unit sphere and 

£ = V^V, V = l, V:=(d/d&,...d/d( d ). 
si 

This is the simplest non-trivial example of flow in a non-homogeneous medium in higher dimensions 
that is related to QAIS. Similar results can be obtained for more complex cases of QAIS, such as 
those related to rational or soliton solutions of KDV hierarchy or finite-reflection groups. 

Proposition 1 Let <p be any solution of VC;f 2 V0 = regular in the closure of the ball 

n : (Ci - CO 2 + (C2 - Q 2 + ■ ■ ■ + (Cd - O 2 < r 2 . 

Then 

where is volume of Q, and 

where integration is taken over surface of d — 1 dimensional sphere dQ centered at £ = C' , and 
s d _i denotes area of this sphere. 

This is an analog of mean value theorem for harmonic functions. 

Proof: is by direct solution of the free-boundary problem (j3J) in d dimensions. Below we 
present solution for all dimensions except d = 2 and d — 4. The proofs for the above dimensions 
differ by technical details only. 

Let us construct a solution assuming that operator q in (124"]) is of the first order. We may write 
down solution of (I24p in the following form, 

P(C) = aG((, + + *(C; CO (25) 

where a and (3 are (^-independent parameters, CO i s a homogenous solution, i.e. C[^/((; (')] = 
0, and G is a fundamental solution 

C[G((, CO] = -&)■■■ S(Q - O- (26) 
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Let us also introduce fundamental solution of the Laplace operator in d- dimensions 

A[G (C,O] = *d-i*(Ci-Ci) • • -%-Cd), Go = p 2 ^, p = ^(Ci - CO 2 + ... + (&- C) 2 (27) 

The fundamental solution G of £ can be found applying Hadamard's expansion, following proce- 
dure developed for QAIS in [3J [5]. Here we briefly describe it: 
Consider a QAIS with the Hamiltonian 

H = A + 2A[logr(C)] = £(C)V£(Cr 2 V£(C) 

and the corresponding free-boundary problem is governed by the elliptic operator 

c = V£(C)- 2 V, (28) 

where r(£) is a polynomial in Ci, • • • , Cd- According to [U |5], the fundamental solution fl26|) can be 
written explicitly as follows 



n(t t>\ — £(0£(C) V^-r 

G K,C) - — ^7^2^ 7i 



i=0 



A- l [G (C,CO] 



(29) 



where n = deg(r(C)), and operators 71, acting on A *[Go], are defined recursively as 

T i+1 = T i A-HT i , r = r(C). (30) 

When % vanish for i > n, the Hadamard's expansion fl29l) truncates, which is a property of QAIS. 
Note that, as follows from the last equation, T n is nothing but the intertwinning operator and 
one can also consider recursive procedure (I3"U|) as a test for algebraic integrability of a system 
with given r(Q (or £(£)) accompanied by simultaneous construction of the intertwinning identity 
HT n = TnA. 

In our case r(() — = d, n = 1, and 

% = Ci, 7i = 2— — . 

d(i Ci 

Taking into account that in d dimensions (except d = 2 and d — 4) 

A" 1 [G ](p) = ^ i yG (p), 



we get 



Substituting 



a: 



G = ( CiCi - ^ ) P 2 - d - (31) 
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into (125]) . and taking (I3TT) into account, we get 

It follows that P vanishes at the boundary p = |C ~ CI = r an d therefore P is a solution of 
free-boundary problem for expanding sphere. From (|25l) . ( 126]) we get 

mio = (i - - co . . .5(0 - oi- 

Applying procedure of derivation of quadrature identities, given in section [3] and taking the above 
equation into account, we complete the proof of the Proposition. 

Similar propositions could be proved for other <i-dimensional analogs of planar integrable el- 
liptic growth systems mentioned in this paper. 

In the planar case, rather than using conformal mappings, one could define the polynomial 
domains as the Laplacian quadrature domains formed by fluid injected through a multipole source 
into an initially empty homogeneous medium. Such a definition of the polynomial domains, as 
Laplacian domains, is acceptable in any dimension. 

It is, however, unclear if higher- dimensional analogs of theorems considered in this paper are 
also valid for the d > 2 Laplacian quadrature domains other than balls. 

In general, proofs use the following procedure: By Green's theorem 

f [p(oc[G](c,o-G(c,oc[p](o)d d c= I (p(oacr 2 dnG(c,o-G(c,oacr 2 9nP(o)ds 

and by (j2J), ©, ([24]), (J2E& we have 

P(C) = OT(C,0- / G(c,CMCK(C)^. 

Jan 

Taking again into account the boundary condition Q, we get the equation for (constant) coeffi- 
cients of q: 

q[G]((,C)= [ G(C,CMCK(C)rf5, CCeda (33) 
Jan 

To further simplify solution of this system one has to use (129]) and properties of Laplacian domains 
(e.g. v n = dr(t)/dt in the case of the sphere). 

6 Conclusion and open questions 

In this article we have considered new applications of the theory of Quantum Algebraically Inte- 
grable Systems (QAIS) to the fluid dynamics of free boundary viscous flows. In two dimensions, we 
used conformal mapping techniques in combination with conservation laws in the form of quadra- 
ture identities to solve whole classes of inverse free-boundary problems. In higher- dimensions 
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we applied a "brute force" approach solving free-boundary problems directly. This approach re- 
lies on methods developed in the theory of QAIS, which limits us to a quite restricted class of 
problems. It will be interesting to find more general and perhaps less technical approaches for 
higher-dimensional cases. 

It is worthy of mention related classification problems, namely the problem of description of 
all elliptic operators (or media, in the language of fluid dynamics) for which polynomial evolution 
of the free-boundary ([9]) can be sustained by a finite number of sources. This question seems to 
be connected with an open problem of classification of QAIS [3l [7J |6]. 
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